Submitted to PRB 



Wavelength Tunability of Ion-bombardment Induced Ripples on 

Sapphire 

Hua Zhou, Yiping Wang, Lan Zhou, and Randall L. Headriclf 

Department of Physics, University of Vermont, Burlington, VT 05405 USA 

Ahmet S. Ozcan, Yiyi Wang, Gozde Ozaydin, and Karl F. Ludwig Jr. 
Department of Physics, Boston University, MA 02215 USA 

D. Peter Siddons 
National Synchrotron Light Source, 
Brookhaven National Laboratory, Upton, NY 11973 USA 
(Dated: February 6, 2008) 

Abstract 

A study of ripple formation on sapphire surfaces by 300 - 2000 eV Ar + ion bombardment is 
presented. Surface characterization by in-situ synchrotron grazing incidence small angle x-ray 
scattering and ex-situ atomic force microscopy is performed in order to study the wavelength of 
ripples formed on sapphire (0001) surfaces. We find that the wavelength can be varied over a 
remarkably wide range - nearly two orders of magnitude - by changing the ion incidence angle. 
Within the linear theory regime, the ion induced viscous flow smoothing mechanism explains the 
general trends of the ripple wavelength at low temperature and incidence angles larger than 30°. In 
this model, relaxation is confined to a few-nm thick damaged surface layer. The behavior at high 
temperature suggests relaxation by surface diffusion. However, strong smoothing is inferred from 
the observed ripple wavelength near normal incidence, which is not consistent with either surface 
diffusion or viscous flow relaxation. 

PACS numbers: 68.35.Bs, 61.10.Eq, 81.16.Rf, 81.65. Cf 
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I. INTRODUCTION 



Energetic particle bombardment on surfaces is known to produce one- dimensional (rip- 
ples or wires) and zero- dimensional (dot) structures at the submicron or nano-scale by a 
self-organization process. Recently, significant experimental and theoretical effort has been 
expended in order to develop ion bombardment patterning methods for the production of 
periodic nanostructures on various substrates.— i 2 i 3 i 4 i 5 i 6 i 7 i 8 These studies have demonstrated 
the potential to tailor surface morphology and related surface properties for novel optoelec- 
tronic and spintronic applications.^*^ In addition, recent work has provided new insight into 
the mechanisms of the instability-driven self-organization process.— 

A significant milestone in our understanding of the origins of a self-organized ripple 
topography formed by ion sputtering is the work of Bradley and Harper (BH) in which they 
proposed a linear continuum equation to describe the main features of ripple formation.— 
The main idea of BH is that smoothing and roughening processes have different wavelength 
dependence, leading to a preferred wavelength where the surface amplitude grows the most 
rapidly. 

However, certain experimentally observed features, such as the saturation of the ripple 
amplitude and the appearance of kinetic roughening are not predicted by the linear 
BH theor y 13 ! 14 An extension of the linear BH theory into the non-linear regime has 
been proposed in order to avoid these shortcomings,— resulting in a noisy version of the 
Kuramoto-Sivashinsky equation: 



The Vo term, which represents the average erosion rate of the unperturbed planar surface can 
be neglected in Eq. (1), since it does not affect the process of ripple formation. The surface 
height h is then in a coordinate system that moves with the average surface during the erosion 
process. r)(x,y,t) is a Poisson noise term related to random fluctuations, uncorrelated in 
space and time, in the flux of the incoming ions. 

Within this theory, ion sputtering produces periodic modulated features (correlated lat- 
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eral ordering) that arises from a competition of a roughening instability mechanism and 
surface relaxation. A roughening mechanism that often dominates the surface morphology 
is curvature-dependent sputtering, which is based on the linear cascade approximation first 
proposed by Sigmund.— However, certain compounds, such as GaSb and InP,^*i& and ele- 
mental materials with refractory-metal seeding,— may also exhibit island agglomeration of 
excess elements by a process related to preferential sputtering. 

In contrast, a wider range of relaxation mechanisms have been proposed in order to ex- 
plain various experimental observations: (i) Surface diffusion (SD) mediated smoothing has 
been proposed to explain the temperature and ion flux dependence of ripple wavelength in 
the high temperature regime.^*^ (ii) The surface erosion smoothing (SES) model in which 
higher-order terms of the erosion process produce a smoothing which mimics surface diffu- 
sion. 8 (iii) The ion-enhanced viscous flow relaxation (IVF) model, which considers surface- 
confined viscous flow driven by surface tension as the dominant smoothing mechanism for 
any material with a disordered surface layer with reduced viscosity— Mayr et al. have em- 
phasized the role of flow of point defects and of thermal-spike induced local melting in the 
mechanisms of radiation- induced viscous flow in the 0.1-1 KeV range. 20 However, determin- 
ing which relaxation mechanism (SD, SES, IVF, etc.) is dominant for various amorphous 
or crystalline substrates, surface temperatures, and ion beam parameters, still remains the 
subject of vigorous research. 

II. LINEAR THEORIES 

In the theoretical approach of Makeev and Barabasi, 11 a general continuum equation 
which describes the evolution of surface morphology during ion sputtering is proposed. Eq. 
1 incorporates the major features of ion-bombardment induced ripple formation and kinetic 
roughening. In the early stage of ion sputtering, Eq. (1) reduces to a BH-type linear theory 
when X x =X y =0, £ x = ^=0. In this article, our discussions of the wavelength tunability of 
ion-bombardment induced sapphire ripples are emphasized within the linear regime. The 
linear terms with coefficients u x and u y represent the curvature dependent ion erosion rates, 
and K xx and K yy are coefficients representing the surface smoothing terms. 

Linear stability analysis indicates that the establishment of a periodic ripple structure 
across the surface depends on the balance between the curvature dependent roughening 
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and surface smoothing mechanisms. 12 Two modes of rippled morphology can be induced 
by ion bombardment, with ripple wavevectors parallel or perpendicular to the projection 
of the ion beams. Regardless of the respective smoothing mechanisms (SD, SES, IVF), 
the wavelength of ion sputtered ripples with orthogonal orientations i x (parallel) and i y 
(perpendicular) are generally expressed as: 



The minimum between £ x and £ y determines which orientation dominates the ion-induced 
ripple topography. In Eq. (2) and (3), the coefficients u x , v y for curvature dependent 
roughening terms are given (following Ref. 11) by : 




(2) 



(3) 



v x = Fd^ {2dts* - duyc 2 + d 2 jy c 2 - dy} , 



(4) 



v v = —Fd 




(5) 
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In the expressions above, the terms are defined as: 




s = sin(0), c = cos(0), 



F = 



f = dy + dy, 

Je P d e -d%dP^/2f 



(6) 



where F is a coefficient relating to the local sputter yield Y{6) as: 



F 



JY(9) 
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In addition, d is the ion energy deposition depth, a and \x are ion energy distribution widths 
parallel and perpendicular to the incoming ion beams, J is the ion flux per area, p is a material 
constant depending on the surface binding energy Uo and scattering cross-section,^ and n 
is the atomic density of the substrate. 

For the SD smoothing mechanism, thermally activated surface diffusion induces surface 
smoothing during ion sputtering. If surface self-diffusion is isotropic, then 

Kxx = K yy = K S d = nf 1 ^ ■ (8) 

Here, D s is the surface self-diffusivity, which has an Arrhenius temperature dependence 
D s =D exp(— AE/k B T). The surface tension (surface free energy per unit area) is 7, and p 
is the areal density of diffusing atoms.— 

For the SES model, the erosion smoothing process is assumed to dominate over 
other smoothing mechanisms. Since it is not temperature dependent, it is thought to 
dominate at low temperatures in some cases. Here, K xx =K XXjSES and K yy =Kyy iSES , 
which are anisotropic with respect to the direction of the incoming ion beams. The coeffi- 
cients representing surface erosion smoothing K xx S es an d K yy ,sES are given (in Ref. 11) by : 



K xx ,ses = {-ml* 2 } + 4s')f + 4c 2 (3/ 2 + 6dyf + d^)f} 



+ F Y4jB W - W 15 W + 10^V + ^V)} , (9) 



d 3 c 2 3dl , . 

Kyy,SES = ^J^f ■ (™) 



For the IVF model, the ion-enhanced surface viscous flow within a thin ion-damaged layer 
dominates the surface smoothing. Hence, 

K XX = Kyy = Kj VF = . (11) 

Vs 
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Here, the ion-enhanced surface viscosity 7] s and the surface tension are assumed to be con- 
stant and isotropic. The depth of the damaged layer is taken to be equal to the ion pene- 
tration depth dM 

In the following sections, the data analysis and discussion of the wavelength tunability of 
ion-bombardment induced sapphire ripples are based on Eqs. 2-11 for the ripple wavelength. 
Additional smoothing mechanisms beyond SD, SES and IVF are considered in section V. 

III. EXPERIMENTAL 

The ripples are produced on sapphire (0001) in a custom surface x-ray ultra-high vac- 
uum chamber (10 -9 torr base pressure) installed at x-ray beamline X21A1 at the National 
Synchrotron Light Source (NSLS) at Brookhaven National Laboratory. More details about 
this real-time surface x-ray characterization facility for dynamic processing are described 
elsewhere.— Ion beam sputtering is performed by either a 3-grid RF plasma ion source or 
a Phi model 04-192 sputter ion gun. The RF plasma ion source is operated at ion energies 
ranging from 300 eV to 1000 eV at a flux of 2xl0 14 ions/cm 2 /sec with a background Ar + 
pressure 4xl0 -4 torr. The irradiation size of the ion beam with a uniform flux (3 cm) is 
large enough to cover the entire sample surface. The Phi sputter ion gun is operated at ion 
energies from 500 eV to 2000 eV at a flux of lxlO 13 ions/cm 2 /sec with a background Ar + 
pressure lxl0~ 4 torr. The sample surface temperature is adjusted from 300 K to 1050 K, 
and monitored by an infrared pyrometer. The chamber is also equipped with reflection high 
energy electron diffraction (RHEED), which is used to determine surface crystallinity. A 
contact mode Digital Instruments Nanoscope-E AFM is used for ex-situ surface morphology 
imaging. 

The x-ray flux after the Si (111) monochromator crystal is 2xl0 12 photons/sec at a wave- 
length of A=l. 192 A with a beam-size of 0.5 mm x 0.5 mm. In the schematic representation 
of the x-ray measurement geometry shown in Fig. 1, the z axis is always taken to be normal 
to the sample surface, and the y axis along the projection of the incident x-ray beams onto 
the surface, ki and kf are the wave vectors of the incoming and scattered x rays, respectively. 
The components of the scattering momentum transfer Q = kf — ki can be expressed by the 
glancing angles of incidence (a^) and exit (a/) with respect to the surface (x-y plane), and 
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the in-plane angle ip. 

In the GISAXS geometry, the incident or exit x-ray beams are fixed near the critical 
angle for total external reflection (0.2° for sapphire). A 320-pixel linear position sensitive 
detector (PSD) is positioned along x axis at the angle ctj with respect to the surface, in 
order to collect in-plane scattered x rays. In terms of scattering momentum transfer, the PSD 
acquires a range of Q x at a constant Q z and Q y (typically, Q 2 =0.92 nm^ 1 , AQ x ^il nm -1 
and (3j/=3xlCr 3 nm" 1 ). Time-resolved GISAXS provides access to the evolving wavelength, 
shape and amplitude of surface ripples. For a full reciprocal mapping (Q x vs. Q z ), an «j = af 
reflection mode scan is performed. 

IV. RESULTS 

In this part, subsections A through D describe a series of systematic investigations of 
the dependance of sapphire ripple characteristics (wavelength, orientation), on experimental 
parameters, including ion energy, ion incidence angle and temperature. 

A. Ripple Formation by Ion Sputtering 

Fig. 2 displays a Q x vs. Q 2 GISAXS scan of the rippled sapphire surface after 45° 
off- normal 500 eV Ar + bombardment. First and second order satellite streaks are readily 
observed, indicating a periodic surface morphology along the x-direction. In this image, 
the specular reflected x-rays are cut off by a narrow beam stop attached to the PSD. This 
provides better sensitivity to off-specular diffuse scattering pattern. 

Fig. 3 shows an AFM image of the as-irradiated sapphire surface. A ripple morphology 
is clearly visible. A wavelength of 32 nm is obtained from the image, which corresponds well 
with that acquired from GISAXS in Fig. 2. 

Fig. 4 shows real-time GISAXS intensities plotted versus scattering momentum transfer 
Q x . The scans are equivalent to a linear slice of the image in Fig. 2 at a constant vertical 
component of the scattering momentum transfer Q 2 =0.92 nm -1 , except that that an offset 
condition («j 7^ «/) was used so that Q y =3x 10~ 3 nm -1 . This was done in order to eliminate 
the need for a specular beam stop. Scans are shown at 10-minute intervals during 45° off- 
normal 500 eV Ar + bombardment. At time t=0, the initial roughness of the sapphire 
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surface only produces a single peak in the diffuse scattering (circles). Two satellite peaks, 
resulting from lateral correlated roughness become visible after 10 minutes. The ripple 
wavelength, i x = 2tt/AQ x , remains almost constant during ion irradiation, but the peak 
intensity continues to increase as a result of an increase in ripple amplitude. 

It is also clear in Fig. 4 that the two satellite peaks develop in an unequal way as 
irradiation proceeds. After 30 minutes, the peak on the positive Q x side is noticeably larger 
than the one on the negative Q x side. At 40 minutes, the larger peak is several times more 
intense than the smaller one. This diffuse intensity asymmetry was also observed in the 
GISAXS study of ion-eroded SiC>2 by Umbach et al.— 

Fig. 5 shows an asymmetric saw-tooth model profile, which is used as a simplified approx- 
imation to the ripple shape. Here, the parallel component of the incident ion beam is along 
the -x direction, as defined above. We note that the term dh/dx in Eq. (1), representing 
the surface local slope, has opposite signs at different sides of the solid saw-tooth. Thus, 
an off-normal incidence will produce different erosion rates on positive and negative slopes. 
The model predicts that the unbalanced erosion will make ripples move like waves across the 
surface in a direction opposite to the projection of the incident beam along the surface.— 
However, we note that a recent study of ripples formed on ion-bombarded glass surfaces 
showed forward propagation of ripples.™ Nonlinear terms of the form (dh/dx)(d 2 h/dx 2 ) 
produce the asymmetric shape, which is observed in our x-ray diffuse scattering measure- 
ments.— Therefore, the appearance of an asymmetric GISAXS pattern indicates the onset 
of this lowest order non-linear term. 

B. Ripple Wavelength Variation with Ion Energy 

Fig. 6 shows the observed dependence of ripple wavelength £ x on ion energy e for ion 
sputtered sapphire at low temperature 300 K (a) and high temperature 1000 K (b), respec- 
tively. This series of sapphire ripples are obtained at 45° off-normal ion incidence. In Fig. 
6(a), square and circle symbols represent the wavelength of ripples produced by the high- flux 
RF plasma ion source and the low-flux ion gun, respectively. The sapphire ripple wavelength 
increases with ion energy at low temperature, which is consistent with observations for ion 
eroded SiO"2, GaAs and Si surfaces . 1 ?i 24 i 2 ^ The data obtained from high/low flux ion sources 
overlap within experimental error at both 500 eV and 1000 eV, which indicates that the 
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ripple wavelength is independent of incident ion flux at low temperature. A non-linear least 
squares fit gives a power law coefficient with p=0.71 for the dependence of the wavelength on 
ion energy (£ x ~ e p ). Also plotted are curves corresponding to p=l and 0.5 for comparison. 
In Fig. 6(b), the ripple wavelength decreases with ion energy. A non-linear least squares fit 
gives a power law coefficient with p=-0.44. Also plotted are curves corresponding to p=-0.25 
and -0.75 for comparison. 

A general formula for a low temperature ripple wavelength along the dominant orientation 
x-axis can be expressed from Eq. (2), (4) and (9), based on the SES model (section II). 
Taking \u x \ ~ Fd and K XXj ses ~ Fd 3 , we can obtain the dependence of the wavelength on 
ion energy as: 



The dependence of d on e is quantitatively accessible with the aid of the ion-collison simulator 
SRIM.^ It indicates that d varies as e v with p=0.48 for a-sapphire at the incidence angle 
of 45°. The p=0.48 obtained from Eq. (12) matches the observed wavelength-ion energy 
relation in Fig. 6(a) reasonably well. However, a quantitative analysis (as detailed in section 
IV-C) by Eq. (12) predicts values of the ripple wavelengths that are an order of magnitude 
smaller than our measured values of £ x , indicating that the SES, which contains no adjustable 
parameters, does not account for the observed ripple wavelength at low temperature.— 

A specific expression for the ripple wavelength based on the ion-enhanced surface viscous 
flow (IVF) model,-^ can be derived from Eq. (2), (4), (7) and (11) as an extension to Eq. 
(12). Inserting Kjvf = ^d 3 /r] s , \u x \ ~ Fd and F=JY(#)/nc into Eq. (2), we can get the ion 
energy dependence for the IVF model as: 



where the coefficient K xx of the IVF smoothing term replaces K xx in Eq. (2). Simulations 
from SRIM at 45° incidence indicate that Y(9) varies as e 5 with 5=0.62. It is assumed from 
previous work on low energy ion irradiation^ that the dependence of ion-enhanced surface 
viscosity on ion energy takes the form Jr] s ~ e~ a with a ~ 1. Hence, l x ~ e p where p=0.67 




(12) 




(13) 
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for sapphire ripples. Here, the sapphire atomic density n is 1.17xl0 29 atoms/m 3 , the surface 
tension of sapphire 7 is 0.91 J/m 2 (Ref. 28), and the ion flux is 2xl0 18 ions/m 2 /sec. The 
ion-enhanced surface viscosity i] s is estimated to be 7.5 xlO 9 Pa-s at an ion energy of 600 
eV. This value is chosen in order to give quantitative agreement between the model and the 
experimental values of l x . 

The fitted power law coefficient of p=0.71 in Fig. 6(a) and that of the IVF model (p=0.67) 
for the dependence of the wavelength on ion energy at low temperature are indistinguishable 
within experimental error. Moreover, we have observed that the RHEED pattern for the 
sapphire (0001) surface disappears upon ion irradiation at room temperature, confirming 
that ion bombardment amorphizes the surface, or at least induce a layer with a very high 
defect density. Similar behavior of surface amorphization under ion bombardment has been 
noted in the study of ion sputtered Si and InP . 16 ' 29 The key idea of the IVF model is that 
this thin layer can relax by a collective motion ("flow"), driven by surface tension. 

The SD model can be useful in predicting the high temperature ripple wavelength. The 
formula for the wavelength with its wavevector along the x-axis can be expressed from Eq. 



where only d and Y(9) are dependent on ion energy. Thus, the SD model gives £ ~ e v where 
p=-0.55 for sapphire ripples at high temperature. The SD model is consistent with the 
energy dependence at high temperature in Fig. 6(b). A more refined model that combines 
both the IVF and SD mechanisms is given in section IV-D, which is also compatible with the 
data in fig. 6(b). Other variations of the SD model that include ion-bombardment effects 
are discussed in section V. We also note that the sapphire surface exhibits a well developed 
RHEED pattern after etching at 1000 K, indicating a higher degree of surface crystallinity 
at this temperature. 

C. Dependence of Ripple Wavelength on Ion Incidence Angle (Low Temperature) 

The observed wavelength-angle phase diagram for sapphire ripples produced by 600 eV 
Ar + bombardment at room temperature is displayed in Fig. 7. The wavelength of sapphire 



(2), (4), (7) and (8), as: 




(14) 
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ripples is varied through a remarkably wide range ( 30 nm to 2 fj,m) by changing the incidence 
angle. Below 40°, the ripple wavelength is particularly sensitive to the incidence angle, while 
the wavelength is relatively constant in the middle range from 40° to 65°. Ion incidence at 
an angle larger than 70° rotates the orientation of the ripples by 90°. 

The theoretical wavelength-angle phase diagram for sapphire ripples produced by 600 
eV Ar + bombardment at room temperature is also shown in Fig. 7. Several variables in 
the expressions of Eqs. (4)-(ll) were calculated using SRIM simulations, including the 
ion penetration depth d, the sputter yield Y(0), and the widths of the deposited energy 
distribution, a and /i. Here, the sapphire surface binding energy Uo is taken to be 2 eV/atom, 
and the ion flux is 2xl0 18 ions/m 2 /sec. We estimate the ion-enhanced surface viscosity to 
be f] s = 7.5 x 10 9 Pa-s, as discussed in Section B. The predicted wavelengths £ x and £ y , 
based on the SES (the dash-dot line and dotted line) and IVF (the solid line and dashed 
line) models, are compared with those experimentally obtained. 

There are three distinct regions in Fig. 7. Region I: v x < v y < 0, K XXj ses > 0, K yy:S ES > 
0, Kivf > 0. This region spans from normal incidence to oblique incidence around 62°. The 
ripple with wavevector parallel to the projection of the ion beam (£ x ) dominates the surface 
morphology in this region. For the SES model, the £ y oriented ripple is predicted to dominate 
the morphology over most of region I, except near the region I/II boundary where the ripple 
orientation rotates by 90°. This is in clear disagreement with the experimental observations. 
On the other hand, the IVF model prediction agrees with the observed ripple orientation 
over the whole range. However, the IVF mechanism does not predict the observed large 
wavelength-angle dependence below 40°. This point will be discussed further in section V. 

Region II: This region is characterized by negative K XXj ses < 0> which prevents the 
appearance of £ x for the SES model, and thus the SES model predicts only £ y ripples in 
this region. In the IVF model, the £ x wavelength increases to infinity at the region II /III 
boundary near 65°. Thus, the IVF model predicts that the dominant ripples will switch their 
wavevector orientation to the y-direction in this region. This boundary could be adjusted 
towards higher angles since it is very sensitive to the change of simulated parameters, such 
as d, d a and d^. The experimental observation is that £ x ripples are still observed, but 
longer scale order in the orthogonal direction begins to build up. Overall, the behavior in 
this region agrees reasonably well with the prediction of the IVF model. 

Region III: v x > 0, v y < 0, K xx>S es < 0, K m> sES > 0, K IVF > 0. The £ x ripple is not 
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stable in either model, since v x becomes an effective smoothing term when it is positive. 
Near 90°, l y either drops to zero (SES) or increases to infinity (IVF). Again, the IVF model 
correctly predicts the observed behavior, at least qualitatively. 

Ex-situ AFM images in Fig. 8 display surface morphologies obtained at different angles 
of incidence for ion sputtered sapphire, corresponding to the observed phase diagram in Fig. 
7. Figs. 8(a) and 8(b) show images for off-normal incidence at 25°, which produces micron- 
scale ripples with wavevector parallel to the ion beam direction, and are readily visible in 
the large-scale image [Fig. 8(b)]. In contrast, 55° incidence, shown in Fig. 8(c), produces 
a well-ordered nanorippled surface with the wavevector parallel to the projection of the 
incoming ion beams along the surface, which has a wavelength similar to that shown in Fig. 
3. We note that at the larger scale in Fig. 8(d) the surface roughness is also correlated with 
wavevector perpendicular to the incoming ion beam, as predicted by l y in the phase diagram. 
Ion incidence at 65° still creates detectable ripples with wavevector parallel to the projection 
of the incoming ion beams in Fig. 8(e), but obvious submicron furrows oriented along the 
ion beam direction are observed in Fig. 8(f). Grazing ion incidence at 75° switches the 
orientation of the ripple wavevector perpendicularly. Fully developed ripples are observed, 
with an unusual rod-like structure, as shown in Figs. 8(g) and 8(h). 

D. Temperature Dependence of Ripple Wavelength 

Fig. 9 shows the observed ripple wavelength dependence on inverse temperature 1/T for 
two different angles of incidence. The ion energy is 600 eV for both angles of incidence. All 
samples are preheated to a chosen temperature and then sputtered at this temperature until 
a well-defined wavelength is established. The ripple wavelength obtained at 45° is constant 
at low temperature and increases significantly when the temperature increase over 700 K. 

We have used the SD mechanism to describes the temperature dependence of the ripple 
wavelength at 45°. The wavelength varies as t ~ (T)~ 1 / 2 exp(—AE/2kBT), where AE is 
the activation energy for surface diffusion and ks is Boltzmann's constant ~ However, this 
expression does not take into account the low-temperature component of ion bombardment 
induced smoothing. Thus in Fig. 9, the observed dependence of the ripple wavelength on 
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temperature at the incidence angle of 45° is modeled with a modified K xx (solid line), as: 

K xx (e,6,T) = K IVF {e,e) + K SD (T) . (15) 

Taking Eq. (15) and the expression for K SD in Eq. (8), we can obtain the temperature 
dependence of the ripple wavelength from Eq. (2). The fitting parameters are the activation 
energy AE for surface diffusion and the areal density of mobile species p involving in surface 
diffusion. An activation energy of AE = 0.72 eV is extracted from Fig. 9 (solid line). Taking 
n=1.17xl0 29 atoms/m 3 , 7=0.91 J/m 2 , and assuming Do with the order of magnitude 10 -8 
m 2 /s, we can estimate p to be of the order of 10 13 m~ 2 . 

The calculated wavelengths for 35° incidence based on Eq. (2), (4), (8), (11) and (15) are 
shown by the dashed line in Fig. 9. A weaker, but still significant temperature dependence 
is predicted, which is not observed experimentally (square symbols). Rather, the experi- 
mental ripple wavelengths at 35° are independent of temperature. This indicates that the 
different smoothing mechanisms (i.e. thermal vs. non-thermal mechanisms) have different 
dependences on the angle of incidence. In particular, the rapid increase in wavelength at low 
angles is inferred to be due to a non-thermal smoothing mechanism that increases at low 
angles to dominate over the other mechanisms, but is not included in our present model. 

Finally, Fig. 10 shows the observed wavelength-angle phase diagram (£ x vs. 9) for 
sapphire ripples obtained at two temperatures, 300 K and at 1000 K. The theoretical curves 
based on Eq. (2), (4), (8), (11) and (15) are also plotted. The solid line for 300 K is identical 
to the solid line in Fig. 7, which is just based on the IVF model. The dashed line for 1000 K 
includes both SD and IVF smoothing effects, combined as shown in Eq. (15). We note that 
the experimental ripple wavelength £ x exhibits a very large increase and is not sensitive to 
thermal activation for incidence angles below 40°. Taken together, these observations again 
indicate a strong non-thermal smoothing mechanism which is not adequately explained by 
any of the models under consideration. 

V. DISCUSSION 

The observations in the section IV-B and IV-C indicate that the IVF and SD models fits 
some of the trends of the wavelength dependence on experimentally accessible parameters. 
However, some of the characteristics during the ripple formation are beyond the current 
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theories as described in sections IV-C and IV-D. Fig. 7 shows that the obtained wavelength 
for incidence angle lower than 40° spans a range of two orders of magnitude from 30 nm to 
2 /im. Figs. 8(c) and 8(d) confirm that ion incidence at 25° only produces surface features 
at the larger scale. Furthermore, ion sputtering at normal incidence does not roughen the 
surface at all from the onset of irradiation, which is confirmed by a real-time GISAXS study. 
This is in contrast to the theoretical behavior which predicts that the ripple wavelength only 
increases slightly as normal incidence is approached. If the nonlinear terms X x and X y are 
introduced into the continuum equation describing the surface motion, kinetic roughening^ 
is expected in the region I of the phase diagram (when \ x ■ \ y > 0), but such roughening is 
not observed. These unusual effects lead us to propose that there is an additional smoothing 
mechanism that dominates the behavior near normal incidence. 

We have considered the fact that the theoretical ripple wavelength is very sensitive to 
the ion range d, so that a small (factor of two) uncertainty in d would have a large (order 
of magnitude) effect on the calculated wavelength for certain models. A factor that is not 
taken into account in SRIM simulations is the ion-channeling effect.^ Ion channeling in 
crystalline surface layers would be strongest near low-index crystallographic directions, and 
could potentially lead to an increase in ion range d near the [0001] axis, and a corresponding 
increase in ripple wavelength. Based on RHEED analysis of sapphire surfaces after ion 
bombardment at different temperatures, there is a striking increase in surface order at 
higher substrate temperatures, and we would expect a stronger channeling effect for the 
experiments at high temperatures than at low temperatures. However, Fig. 10 does not 
show any effect of temperature at the lowest angles of incidence, effectively refuting this 
idea. Furthermore, the cutoff of ion-channeling by a low energy threshold in the range 
between 0.1-1 keV may weaken its possible influence on the formation of erosion patterns.— 

Another important factor that is not taken into consideration by current models to explain 
surface morphology created by ion sputtering is ion impact induced lateral mass redistri- 
bution .^22*22 it can induce a form of surface smoothing process that is different from the 
SD, SES or IVF relaxation mechanisms previously discussed. Impact-induced downhill cur- 
rents have been identified as the driving force underlying the ultra-smoothness of surfaces 
resulting from ion assisted film deposition. 33 The ion induced lateral currents gives rise to 
curvature dependent (d 2 h/dx 2 , d 2 h/dy 2 ) smoothing terms, which can weaken or even cancel 
the curvature dependent roughening term in Eq. (1) [see Eq. (6) in Ref. 4]. No ripples can 
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form under conditions where this type of smoothing term is dominant, since the wavelength 
dependence exactly matches that of the prevailing roughening mechanism. Thus, the idea 
that this additional smoothing mechanism plays a role at low angles of incidence provides 
a route which may explain the anomalous wavelength-angle phase diagram in Figs. 7 and 
10. We also note that the lateral current term is expected to be strongest at low angles of 
incidence.— 

The behavior of the ripple formation at high temperature can be explained reasonably by 
the SD model. However, another important fact observed in our experiments, but not dis- 
cussed above, is that thermal annealing at 1000 K without ion irradiation does not produce 
any distinguishable decay of the amplitude of as-prepared ripples, in contrast with previous 
studies on Si and Ag surfaces . 34 ' 35 This indicates that the surface smoothing at high tem- 
perature is not by a type of surface diffusion that is purely thermally activated. Rather, it 
is likely to involve to the creation of mobile species on the surface during ion bombardment. 
Further work on the flux dependence of the i x at high temperature will assist us in clari- 
fying the dominant creation process for mobile species underlying the ion-enhanced surface 
diffusion.— 

VI. SUMMARY 

In summary, the formation and characteristics of ripple morphologies on sapphire 
surfaces produced by ion sputtering are systematically investigated by in-situ GISAXS 
and ex-situ AFM. The ripple wavelength can be modulated effectively in a wide range of 
20 to 2000 nm by changing the ion incidence angle, ion energy and temperature. This 
phenomenon provides an easy route to fabricate nanostructured surfaces for exploring novel 
nanoscale phenomenon. The IVF and SD smoothing mechanisms are shown to play an 
important role in the formation of the sapphire ripple structure. The possible importance 
of impact-induced lateral currents as a smoothing mechanism should also be investigated 
further. 
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Figure Captions 

FIG I. Schematic of the x-ray scattering geometry. The z-axis is always taken to be 
normal to the sample surface, and the incident x-ray beam defines the y-z plane, ki and kf 
represent the incident and scattered wave vectors, respectively. The components (Q x , Q, y 
and Q 2 ) of the scattering momentum transfer Q = kf— ki are defined by the glancing angles 
ttj, otf and in-plane (x-y plane) angle tp. The linear PSD is oriented with its long direction 
in the plane of the surface. 

FIG 2. Grazing incidence small angle x-ray scattering (GISAXS) patterns acquired from 
Qx-Qz scans. Asymmetric GISAXS patterns indicate the shape anisotropy on the ripples. 
The bar of gray shades represents the logarithmically scaled intensity. 

FIG 3. A 500 nm by 500 nm AFM image taken after the GISAXS measurements. The 
ripple wave vector is parallel to the projected ion beam direction. In the text, this is referred 
to as the £ x orientation of ripples. 

FIG 4. Time-resolved GISAXS measurements indicate the increase of the lateral correla- 
tions on a sapphire surface during ion exposure. The surface ripples are produced by 500 eV 
Ar + bombardment at 45° incidence. The vertical component of the scattering momentum 
transfer Q z was fixed at 0.92 nm" 1 and Q y was offset to 3 x 10~ 3 nm" 1 . 

FIG 5. An asymmetric saw-tooth model profile (a < f3). The parallel component of the 
incident ion beam along the substrate surface is along the -x direction, as defined in Fig. 1. 
The term dh/dx represents the local slope of ion etched surface curvature. The magnitude 
of a and (3 are exaggerated for clarity. 

FIG 6. The dependence of the ripple wavelength (£ x ) of ion etched sapphire on the 
incident ion energy (e). (a) 300 K substrate temperature. The square symbols are for the 
high flux ion source and the circle symbols are for the low flux ion source. The solid line is 
the curve for the best fit exponent of p=0.71. For comparison, curves for p=l and p=0.5 
are plotted as dashed and dotted lines, respectively, (b) 1000 K substrate temperature. The 
solid line is for the best fit exponent of p=-0.44. For comparison curves for p=-0.25 and 
p=-0.75 are plotted as dashed and dotted lines, respectively. 

FIG 7. The wavelength-angle phase diagram for sapphire ripples produced by 600 eV Ar + 
bombardment. Circle symbols: experimental wavelength £ x at room temperature. Square 
symbol: experimental wavelength £ y at room temperature (Only one point at 75° is shown). 
The error bars are comparable to the height of the symbols. Two sets of curves are shown: 
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The upper set is for the IVF model, with the solid line representing the simulated i x wave- 
length, and the dashed line representing the i y wavelength. The lower set of curves is for 
the SES model, where the dot-dashed line represents the l x wavelength and the dotted line 
represents the i y wavelength. Experimentally, the three regions of the phase diagram rep- 
resent the ranges of angles where £ x and £ y are dominant (I and III, respectively), and a 
narrow transition region ( region II). The SES model predicts that £ y ripples dominate in all 
three regions, which is not consistent with the data. The IVF model gives better agreement 
with the data, as described in the text. 

FIG 8. AFM images of surface morphology at different angles of incidence: (a) and (b) 
are for an incidence angle of 25°, (c) and (d) are for an incidence angle of 55°, (e) and (f) 
are for an incidence angle of 65°, and (g) and (h) are for an incidence angle of 75°. Two 
different image sizes are shown for each angle of incidence. The white arrow indicates the 
projection of the ion beam direction along the surface. The progression of ripple orientation 
and wavelength can be seen in this sequence of images. In particular, we note the large scale 
ripples at 25° in (b), followed by much smaller wavelength ripples at 55° in (c) and 65° in 
(e), and finally the rotation of the ripple orientation in (g) and (h). 

FIG 9. The wavelength l x as a function of inverse temperature at different angles of 
incidence and constant ion energy of 600 eV. Theoretical curves based on Eqs. (2), (4), (8), 
(11) and (15) are shown for the set of ripples wavelength i x obtained at incidence angles of 
35° (dash line) and 45° (solid line) off-normal incidence. 

FIG 10. The wavelength-angle phase diagram for sapphire ripples produced by 600 eV 
Ar + bombardment at different temperatures. The error bars are comparable to the height 
of the symbols. The curves for £ x based on Eqs. (2), (4), (8), (11) and (15) are compared 
with the measured ripples wavelengths obtained at 300 K (solid line) and 1000 K (dash line), 
respectively. 
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